Stable Bloch oscillations of cold atoms with time-dependent interaction 
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We investigate Bloch oscillations of interacting cold atoms in a mean-field framework. In general, 
atom-atom interaction causes dephasing and destroys Bloch oscillations. Here, we show that Bloch 
oscillations are persistent if the interaction is modulated harmonically with suitable frequency and 
phase. For other modulations, Bloch oscillations are rapidly damped. We explain this behavior in 
terms of collective coordinates whose Hamiltonian dynamics permits to predict a whole family of 
stable solutions. In order to describe also the unstable cases, we carry out a stability analysis for 
Bogoliubov excitations. Using Floquet theory, we are able to predict the unstable modes as well as 
their growth rate, found to be in excellent agreement with numerical simulations. 
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The dynamics of quantum particles in periodic poten- 
tials subjected to uniform acceleration has a long his- 
tory and fascinating physics. Electrons in tilted peri- 
odic potentials present dynamical localization and may 
undergo coherent oscillations, both in real and in mo- 
mentum space, known as Bloch oscillations (BOs) 
BOs were observed for the first time as coherent oscilla- 
tions of electronic wave-packets in semiconductor super- 
lattices @, IH • These oscillations are due to interference of 
partially scattered wave amplitudes and therefore can be 
observed for any phase-coherent waves accelerated in pe- 
riodic potentials. BOs have been directly observed with 
ultracold atoms 0, 0| , Bose-Einstein condensates [f| and 
noninteracting fermions in tilted optical lattices. 

After excitation BOs persist until the quantum parti- 
cles lose their phase coherence. In Bose-Einstein conden- 
sates interatomic interactions lead to a rapid broadening 
of the momentum distribution and a strong dephasing. 
Even in the most favorable experimental conditions, only 
few cycles are usually observed. It is then believed that 
nonlinearities must generally lead to a breakdown of the 
BOs d SOU. 

In cold-atom experiments, one is able to change the 
interaction strength by means of Feshbach resonances 
111 Il2fl. Reducing the interaction to zero, Gustavsson 
et al. [13| were able to increase the dephasing time of 
a cloud of Cs atoms from a value slightly larger than a 
Bloch period Tb to more than 2 x 10 4 T-q. A different way 
of obtaining stable Bloch oscillations was proposed by 
Salerno et al. LJ] : properly designing the spatial depen- 



dence of the scattering length around the zero crossing 
in a Feshbach resonance should result in long-living BO 
of bright solitons. In this Letter, we propose instead to 
modulate the interaction harmonically in time by an os- 
cillating, but spatially homogeneous magnetic field close 
to a Feshbach resonance [l5[ , which should be much eas- 
ier experimentally. Although in general atom-atom in- 
teractions lead to a rapid dephasing of BOs, remarkably 
with some harmonic modulations we find stable oscilla- 




FIG. 1: Condensate density |^ n (t)| 2 for two interactions 
modulated harmonically at the Bloch frequency lob = F: 
(a) g(t) = g cos(Ft) and (b) g(t) = g sin(Ft) with F = 0.2 
and go = 1. Data from numerical integration of the Gross- 
Pitaevskii equation JT}. The initial width of the wave packet 
is (jo = 10, equal to the amplitude xb ~ 2/F of free BOs. 
While the BOs in (a) are perfectly stable, they are rapidly 
destroyed in (b). 



tions of the condensate [see Fig. [TJa)], while for others 
the oscillations are rapidly dephased [see Fig. [Tfb)] . 

Our starting point for this study is the mean-field 
Gross-Pitaevskii equation for the discrete tight-binding 
model [l6| 

= - <F ra -i + Fn% n + .g(t)|*„| 2 *„ , (1) 

where ^ n denotes the wave function at each potential 
well n of the optical lattice. The overdot indicates the 
time derivative, the hopping parameter is set as the unit 
of energy, and h = 1. We assume that initially the atomic 
cloud has a Gaussian density profile centered at rest in 
the lattice, * n (0) = (27rcr 2 )" 1/4 exp [-n 2 /4cr 2 ]. In the 
noninteracting case git) = 0, the atomic cloud oscillates 
with the Bloch frequency lob — 2it/Tb = F. 

To investigate the interacting case, we first solve ([1]) 
numerically by means of the fourth-order Runge-Kutta 
method. Figure Q] shows the time evolution of the con- 
densate density j^n^)! 2 for two different harmonic inter- 
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actions modulated at the Bloch frequency. For g(t) = 
go cos(Ft), we observe that the initial Gaussian shape is 
preserved over time, and the condensate performs per- 
fectly stable BOs with frequency u>b- On the contrary, 
for g(t) — g sm(Ft), the initial shape is distorted al- 
ready after few Bloch cycles. As a consequence, BOs are 
rapidly damped. 

In order to understand how the interaction affects BOs 
of a Gaussian wave packet, we start with a descrip- 
tion in terms of collective coordinates ||. As free vari- 
ables, we chose the center-of-mass position x(t) = (n) = 
E„ n l*n(*)| 2 and width w(t) = {(n - x(t)) 2 ) and their 
conjugate momenta p(t) and b(t), respectively. These are 
defined by their generating role, —id p ^ n = n^ n and sim- 
ilarly for b. The wave function is up to a global phase 
parametrized as 



y n (t) = e^ n A(n-x(t) lW (t),b(t)), 



(2) 



where the rapidly varying exponential exp[ip(t)n] is fac- 
torized from a smooth Gaussian envelope 



A(n, w, b) 
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The equation of motion (JTJ) derives as i^f r 
from the nonlinear Hamiltonian 



(3) 
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Inserting the ansatz (J2]), Taylor-expanding the discrete 
gradient to second order, and performing the Gaussian 
integration, we find the effective Hamiltonian 



H r ,c = FX 
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By construction, the collective coordinates obey the 
canonical equations of motion 
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(6a) 
(6b) 

(6c) 
(6d) 



The initial conditions for the case under study are 
x(0) = 0, p(0) = 0, w(0) = erg, 6(0) = 0. Equation 
([6]) shows that within this ansatz, the momentum of the 
centroid is always given as p(t) = —Ft which then serves 
as the driving term in the other equations. The au- 
tonomous equations (16c]) and (|6dj) describing the width 
can be solved first, their solution finally conditioning the 
centroid motion obtained by integrating (|6b[) . 



Let us study the implications of these equations in 
some simple cases. The complex width of the Gaussian 
envelope, z = [w^ 1 — 4i6] _1 , evolves according to 



z = i cos Ft 
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In the linear case g — 0, this equation has the exact 
solution Rezo == w(Q) — a 2 , and lmzo(t) = F^ 1 sin Ft. It 
describes a Gaussian wave packet with a breathing width 
w a {t) = a 2 [l + (sin Ft) 2 /(Fa 2 ) 2 ] . This solution is valid 
for small 1/(Fctq), i.e., sufficiently broad wave packets. 
Within the collective-variable approach, the amplitude 
of BOs is actually determined by R :— (8Rez) _1 = (1 + 
16b 2 w 2 )/(8w), appearing in (|6b)) . In the linear case, its 
constant value i?o = o"q^ 2 /8 <C 1 yields the usual BOs 
xo(t) — xb cos Ft with amplitude xb = 2(1 — Rq)/F. 

A constant nonlinearity g(t) — go leads to a damping of 
the oscillation amplitude. With the highly nonlinear and 
rapidly oscillating equations of motion (|6c| and (|6d[) one 
finds R = g W(t) cos Ft + g 2 /(l6Trw(t) 2 ). The first term 
contains contributions from w(t) and b(t) that are very 
effectively suppressed over a Bloch cycle by the oscillating 
factor cos Ft. The decay to lowest order in go therefore is 
driven by the second term, and the amplitude decreases 
initially like x B (t) = (2/F)(l - R - 2g%R 2 l t 2 /Tr) @,[IrJ. 

Let us now consider a time-dependent interaction of 
the form g(t) = go sm(ujt + 5). The equations of motion 
for the collective coordinates allow to identify a whole 
family of values (lo, 5) that result in strictly periodic mo- 
tion. Formally, we can combine (|6c[) and (|6d[) to a single 
equation v = X(v,t) for v = (w,b). This Hamiltonian 
flow is driven by a vector field X(-,i) that depends ex- 
plicitly on cos Ft and g(t). Let now u> = (y\lv-i)iOB with 
V\ , f2 G N be commensurate with the Bloch frequency 
lu b = F. Then X(-,t + T) = X(-,t) is periodic with 
T = v%Tb- Since different trajectories v(t) cannot inter- 
sect, the flow dynamics is strictly periodic if there exists 
a time r such that v(t - T/2) = v(t + T/2). This in 
turn is guaranteed by the equation of motion if X(-, t) is 
odd around r, X(-,t — t) = —X(-,t + t). But since X 
contains cos Ft, this time can only be one of its zeroes 
Tj = (2j + l jjF^ 1 , j E N. This then requires also the in- 
teraction to be odd at that point, which fixes Sj = —uTj. 
Therefore, the family of time-dependent interactions with 
periodic solutions contains all linear combinations of 



g(t)=g sm I -±[Ft- -(2j + 1) 



v 2 



(8) 

for example cos ((2n + l)Ft), s'm(2nFt), cos(Ft/3) or 
s'm(2Ft/3). In all such cases, the periodicity of w(t) 
and bit) predicts, via (|6b[) . perfect BOs in spite of the 
underlying nonlinearity. 

This stabilization mechanism does not rely on sup- 
pressing the nonlinearity at the band edges, which can in 
general be helpful to avoid Landau-Zener tunneling [l6| 
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FIG. 2: Centroid x(t) for (a) g(t] = 30 cos(wi), (b) g{t) = 
go sin(tdt) with lo = F — 0.2 and go = 1.0. The solid line 
is obtained by numerical integration of Q; circles connected 
by dashed lines are the solution of the collective-coordinate 
Eqs. (|6a|) - (|6cJ) ; In the stable case (a) both curves agree, 
whereas in (b) the collective coordinates predict a later onset 
of damping. 



but is irrelevant in our single-band description. Instead, 
the simplest stable interaction of type g(t) = ± cos(Ft) = 
±cos(p(i)) has its maximum amplitude with positive 
(negative) sign at the band center p(t) = and nega- 
tive (positive) sign at the band edge p(t) = n. Clearly, it 
does not vanish at the band edge, but rather at the band 
midpoint, where the mass changes sign and dynamical 
instabilities would start to appear 0, [ItJ ■ 

We have tested the collective-variable predictions by 
comparing the centroid x(t) as obtained from equa- 
tion ([6]) with the results from the numerical integration 
of fTJ. Collective coordinates provide an excellent de- 
scription of long- living BOs [Fig. [DJa)]. In the damped 
case g(t) = go s'm(Ft) this approach also predicts a break- 
down of BOs, but not quite correctly [Fig. [U^b)]. This 
ansatz allows only the gaussian width to respond to the 
nonlinear perturbation, but it cannot capture the actual 
dephasing mechanism, which involves the excitation of 
fluctuations as apparent from Fig.[JJ(b). 

To proceed beyond collective coordinates we start with 
an infinitely wide wave packet, namely an atomic cloud 



with uniform density |\& 
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po- In this case, equation 



(HJ) is exactly solved by \l/°(f) = yfpo exp[— i<j)(t)] with 



4>{t) = Fnt - — sin(Ft) 
F 



dt'n{t') 



(9) 



in terms of the chemical potential fi(t) = pog(t). This 
solution predicts that an infinitely narrow momentum 
distribution will perfectly Bloch oscillate in momentum 
space, i.e., the condensate current will oscillate homoge- 



= (N P0 )-ij: k (k-p) 2 \* k f = 

2 ' is given by these ampli- 



neously without signatures in position space, and this for 
arbitrary interaction. 

In the following we study the broadening of the mo- 
mentum distribution due to the growth of small perturba- 
tions. If the wave function deviates slightly from the ho- 
mogenous solution, SSf n (t) = [yfpo + 3>(n, £)] exp[— i<p(t)], 
linearizing {JJ) gives the equation of motion for the fluc- 
tuation Q(n, t): 

i$ = -cos(Fi)$" + 2isin(Fi)$' + 2/i(i)Re$.(10) 

$(n, i) is assumed to be a smooth function with spatial 
derivative = d n $- The interaction term acts on Re<I > , 
i.e., the density fluctuation. Separating real and imagi- 
nary part, <I>(n, t) — s(m, t) + i d{m, t) while transforming 
to the reference frame m = n — Xq (t) of the unperturbed 
BO with Xo(t) = 2F~ X cos Ft, and going to the momen- 
tum representation brings (|10[) to 

4 = - [k 2 cos(Ft) + 2/x(i)] s k , (11a) 
s k = k 2 cos(Ft)d k . (lib) 

These are the Bogoliubov-de Gennes equations for ex- 
citations of a homogeneous condensate with chemical 
potential /i(t) and a time-dependent mass, such that 
e°(£) = k 2 cos(Ft). The width of the momentum- 
space distribution, (Afc) 2 

tudes. Thus, their stability is the key to understand the 
stability of BOs. 

In Eqs. (fTTj) . we have again the structure of a first-order 
Hamiltonian equation of motion (now with a linear evo- 
lution operator) depending on cos(Ft) and /i(t) = pog(t). 
Thus the same argumentation as for the collective coor- 
dinates holds: solutions are periodic and hence pertur- 
bations do not grow if g(t) is of the type {8]). 

But now we are also able to make quantitative pre- 
dictions about both stable and unstable cases. Thanks 
to the linearity of (fTTj) . Floquet theory [HI applies if the 
interaction is modulated at a commensurate frequency 
lo = (v\/v2)u)-r. The stability of this system can be as- 
sessed by integrating the differential equations (fTT|) over 
one period T = v%T^ with two different initial conditions: 
{4(0) = 1,4(0) = 0} and {4(0) = 0,4(0) = 1}. A 
perturbation with wave vector k is stable when |A^| < 1, 
where A k = (1/2) [d\{T) + 4(T)] . In the periodic 
cases Afc = 1. An unstable fluctuation with wave vec- 
tor k grows exponentially with the Lyapunov exponent 

X k =T~ 1 log [A k + v/A^F"! ■ 

Figure [3] indicates by black squares those parameters 
(w/wbi k) in the case g(t) — go sin(wt), for which |Afc| < 1 
assures stability. Full lines correspond to the stable cases 
at lo/ujb = 2n/5, n € N, predicted by (JSJ); other stable 
lines are not resolved due to the finite frequency resolu- 
tion 0.1 used for the plot. For most frequencies, there are 
many fc-vectors belonging to unstable excitations, so the 
BO will decay [Fig.Q] (b)]. 
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FIG. 3: Stability map for g(t) = gosin(ujt), with lo-q — F = 
0.2 and go = 1-0. Black squares indicate the wave vectors k 
of stable excitations at each frequency ui. The stable cases at 
uj/ujb = 2n/5 predicted by ([8| appear as full lines; the other 
stable lines are not resolved due to the finite frequency step 
0.1 used here. For the case of cos(wt), the diagram is similar, 
with stable lines at w/wb = (2n + l)/5. 




FIG. 4: Momentum density |^fc| 2 for selected fc-modes in the 
unstable case g(t) = go sin(Ft). The original wave function is 
centered around k = 0. The solid line marks the growth rate 
of the most unstable mode as predicted by Floquet theory. 
The growth of this mode precedes the damping of the centroid 
motion that sets in at t « 200Tb. Numerical parameters: 
(to = 100, fio = 0.01, F = 0.2. 



For a quantitative test of Floquet-theory predictions, 
we choose a larger width cr = 100 and a weak inter- 
action parameter /i n = 0.01. As obvious from Fig. [3J 
unstable modes are typically encountered at intermedi- 
ate values of k. Floquet analysis predicts as the most 
unstable mode k w 0.44 with a Lyapunov exponent of 
Amax ~ 0.1815 Tg 1 . The numerically calculated centroid 
motion (similarly to Fig. [2]) starts to get damped around 
200Tb. The fc-space picture, Fig. HJ is much more re- 
vealing and shows that the most unstable mode starts to 
grow much earlier. The numerically observed dominant 
mode and also its growth rate are in perfect agreement 
with our Floquet-theory predictions. 

To conclude, we have shown that stable BOs of cold 
atoms are possible when the atom-atom interaction is 



modulated in time with suitable frequency and phase. 
Collective coordinates are found to accurately describe 
undamped BOs. For commensurate frequencies lj and 
lob, we have identified a class of time-dependent interac- 
tions, for which the BOs remain perfectly periodic. But 
collective coordinates cannot capture the dynamics of un- 
stable cases where the main wave packet decays by ra- 
diating excitations. In order to explain the dephasing 
of the oscillations in unstable cases, we applied a stabil- 
ity analysis based on Floquet theory that is in excellent 
agreement with the numerical results. Using harmoni- 
cally modulated interaction thus opens the possibility of 
studying the precise influence of other dephasing mech- 
anisms on Bloch oscillations, such as deliberately intro- 
duced disorder fljl ]. 
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